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Abstract 

We monitored water-vapor masers around the semi-regular variable star 
R Crateris with the Japanese VLBI Network (J-Net) at the 22 GHz band during 
four epochs with intervals of one month. The relative proper motions and Doppler- 
velocity drifts of twelve maser features were measured. Most of them existed for longer 
than 80 days. The 3-D kinematics of the features indicates a bipolar expanding flow. 

The major axis of the asymmetric flow was estimated to be at P.A. = 136°. The 
existence of a bipolar outflow suggests that a Mira variable star had already formed 
a bipolar outflow. The water masers are in a region of apparent minimum radii of 
1.3 X 10^^ m and maximum radii of 2.6 x 10^^ m, between which the expansion ve¬ 
locity ranges from 4.3 to 7.4 km s“^. These values suggest that the water masers are 
radially accelerated, but still gravitationally bound, in the water-maser region. The 
most positive and negative velocity-drifting features were found relatively close to the 
systemic velocity of the star. We found that the blue-shifted features are apparently 
accelerated and the red-shifted apparently decelerated. The acceleration of only the 
blue-shifted features seems to be consistent with that of the expanding flow from the 
star. 
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1. Introduction 


The gas dynamics of circumstellar envelopes (CSEs) during the latest phase of the 
asymptotic giant branch (AGB) stars is one matter of debate to elucidate the physical mech¬ 
anisms of energetic mass loss of evolved stars. The AGB stars evolve in a sequence of Mira 
variables, IRC (or AFGL) objects, OH/IR stars, proto-planetary nebulae (PPNs), and plane¬ 
tary nebulae (PNs) (^00 [cite]cite.TakaTakaba et ah (1994)). Recent HST infrared images (e.g. 
^0 [citejcite.SahaiSahai et ah (1998)) have shown asymmetric or bipolar outflows around the 
central object of PNs; PNs show axisymmetric morphologies in over 50% of the cases (e.g. ^00 
[citejcite.SahaiSahai et ah (1998); ^00 [cite]cite.MancManchado et ah (2000)). On the other 
hand, the CSEs of AGB stars are generally spherically symmetric (000 [cite] cite.MancManchado 
et ah (2000)). The most intriguing question is when the change of morphology occurs between 
the AGB and PN phases. 

Water-vapor masers have been well observed around many evolved stars. A cluster of 
water-maser features exists at 10^^ — 10^^ m from the central star (000 [citejcite.JohJohnston 
et al.(1985); 000 [citejcite.Lane-aLane et al.(1987)). Maser action is predicted in the inner part 
of the CSE, and the size of the maser region increases with the stellar mass-loss rate (000 
[cite]cite.CookCooke, Elitzur (1985)). Distributions of water masers in CSEs are complex, and 
change with time (e.g. 000 [citejcite.LaneLane (1982); 000 [citejcite.JohJohnston et al.(1985)). 
The time variation was reported on a scale ranging from a few weeks to a few years (000 
[cite]cite.Bo-aBowers et al.(1989)). The VLA, VLBA, and MERLIN observations have moni¬ 
tored the water-maser features: the epochs were widely scattered from a half year to several 
years (000 [cite]cite.MarMarvel(1996)). Measurements of the relative proper motion of compact 
maser features with high spatial resolution can be used to investigate the structure of CSEs 
and to trace those asymmetries or bipolarities (030 [cite]cite.MarMarvel(1996)). Measurements 
of 3-D velocities of water-maser features as well as their radial-velocity drifts (000 [cite]cite.Im- 
almai et ah (1997)) enable us to elucidate the dynamics of the mass-loss process in more detail. 
For this purpose, monitoring observations within intervals of a few weeks are necessary to verify 
whether the same maser features are observed within the epochs. 

R Crt is a semi-regular variable star of SRb type in spectral classihcation M7 (000 
[cite]cite.KholKholopov et al.(1987)). R Crt is classihed as a Mira variable from a color-color 
diagram of IRAS observations (000 [cite]cite.TakaTakaba et ah (1994)). The R-band magni¬ 
tude varies between 9.8 and 11.2 mag with a pulsation period of 160 days. ^^CO(J=2—1) and 
(J=l—0) and ^^CO( J=2— 1)(0OO [cite]cite.Kaha-aKahane, Jura (1994)) as well as OH and SiO 
masers (000 [cite] cite. Le-SLe Squeren et al.(1979); 030 [cite] cite. Jew Jewell et al.(1991)) were 
detected at the envelope. The mass-loss rate of R Crt was estimated to be 1.0 x 10“® Mq 
yr“^ (000 [citejcite.Kaha-aKahane, Jura (1994)) from CO observation. The systemic velocity 
for R Crt is 10.8 km s“^ (000 [cite]cite.Bo-bBowers (1992)), which was estimated from the cen- 


2 


Table 1. J-Net telescopes that participated in the observations. 


Station Designation Diameter (m) Ap. efficiency {%) Tsys at zenith (K) 


Mizusawa . M 10 36 140 - 330 

Kashima. O 34 57 280 - 340 

Nobeyama . N 45 63 240 - 290 

Kagoshima. K 6 40 170 - 220 


ter velocity of the velocity coverage of a CO observation profile. Water masers around R Crt 
were first detected by |e]cite.DickDickinson et al.(1973) ([citejcite.DickDickinson et al.(1973)). 
In this paper we present the results of monitoring observations for water masers around R 
Crt. We compared the positions of the water-maser features among epochs separated by a few 
weeks, and measured radial velocities throughout the epochs. Based on the 3-D motions and 
radial-velocity drifts, we discuss the dynamics of the mass-loss flow from this star. 

2. Observations and Data Reduction 

Monitoring observations of R Crt at R.A. (B1950.0) = 10^58™6fl4 and Deck 
(B1950.0) = —18°03'18727 were made by using the Japanese VLBI Network (J-Net, 000 
[cite]cite.OmoOmodaka et ah(1994)) from 1998 March to June, at four epochs with an interval 
of a few weeks. J-Net is comprised of three of the National Astronomical Observatory of Japan’s 
telescopes, including the 45-m telescope at Nobeyama Radio Observatory D the 10-m tele¬ 
scope at Mizusawa, and the 6-m telescope at Kagoshima, as well as one of the Communications 
Research Laboratory’s telescopes, the 34-m telescope at Kashima. Table 0 gives the character¬ 
istics of the telescopes, and table 0 the observations. 

In each epoch R Crt was observed at 5 to 8 intervals for 20 to 40 minutes alternating a 
scan for 5 to 10 minutes toward 3C 273B as a calibrator. The observed signals were recorded 
with the VSOP terminal (000 [cite]cite.KawaKawaguchi et ah(1994)) in one base band channel 
with a band width of 16 MHz in 2 bits per sample. The recorded data were correlated using 
the Mitaka FX Correlator (000 [cite]cite.ChikChikada et al.(1991)). The output of the correla¬ 
tor had an 8 MHz bandwidth with 512 spectral channels, corresponding to a Doppler velocity 
resolution of 0.105 km s“^. The velocity coverage was from Wsr = —16 km s“^ to 36 km s“^ 
centered on the water-maser line. 

We used the NRAO’s AIPS software to calibrate, image and £t the brightness distri¬ 
bution. The strongest water maser component detected at Wsr = 15.03 km s“^ was used as 
phase reference in each set of observations. The synthesized beam of the array was 8 mas x 3 

P) Nobeyama Radio Observatory is a branch of the National Astronomical Observatory, operated by the 
Ministry of Education, Culture, Sports, Science, and Technology!, Japan). 
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Table 2. Time table of the monitoring observation. 


Experiment 

name 

Epoch 

(1998) 

Stations 

Observation 

duration (hr) 

Typical rms thermal 

noise (mJy) 

Number of detected 

features 

j98090 .... 

March 31 

MKNO 

10 

34 

14 

j98116 .... 

April 26 

MKNO 

10 

55 

13 

j98148 .... 

May 28 

MKNO 

9 

70 

12 

j98170 .... 

June 19 

MKNO 

8 

80 

5 


mas with a position angle of — 40° in natural weight. The typical rms thermal noise level in 
the channel maps ranged from 34 mJy to 80 mJy between the hrst and last epochs (see table 
The accuracy of the intensity scale is estimated to be better than 50%. The accuracy of 
the relative angular position ranged from 0.2 mas to 0.9 mas in the R.A. direction and 0.2 mas 
to 1.1 mas in the declination direction, in a single velocity channel. 

A maser spot is a single-velocity component, while a maser feature is a group of maser 
spots within several mas in positions and 0.1 — 0.5 km s“^ in velocities. The dehnitions of 
‘maser spot’ and ‘feature’ used in this paper are the same as those used in |e]cite.lm-blmai 
et al.(2000) ([cite]cite.Im-bimai et al.(2000)) and references therein. These dehnitions are the 
same as those given by |t|e]cite.GwinGwinn (1994) ([citejcite.GwinGwinn (1994)). The positions 
of the brightness peaks in the maser spots were estimated using the AIPS task JMFIT. We 
traced the position and velocity at the brightness peak in each of the features. We determined 
the velocity at the brightness peak of a feature by quadratic htting of the observed intensities 
against the velocities of the strongest three maser spots in the feature. The position of the 
feature at the calculated peak velocity was estimated on the line-connecting positions of two 
velocity-adjacent spots assuming a constant velocity gradient between the two spots. 

3. Results and Discussion 

For proper-motion measurements, the four maps were superposed so as to situate the 
position-reference feature (Wsr = 15.03 km s“^) at the same position in the superposed single 
map. Figure |I] shows the observed proper motions of the maser features in the R.A. and 
Deck directions. We also measured the Doppler-velocity drifts of the maser features that 
are shown in hgure |T[ The uncertainty of the relative velocities of the maser features was 
assumed to be 0.05 km s“^, one half of our velocity resolution. Since the position-reference is 
also moving relatively, the mean proper motion of the features was calculated and subtracted 
from each of the proper motions. Table ^ summarizes the relative positions, proper motions, 
Doppler velocities, Doppler-velocity drifts, and peak intensities of the maser features. A letter 
is assigned to each of the features, and hereafter we use these assignations. We conhrmed that 
the strongest features (D, H, and K) were detected at four epochs, and existed for at least 80 
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900 940 980 


JD-2450000 


Fig. 1. Time variation of R Crt maser features in (a) R.A., (b) DecL, and (c) the Doppler velocity. 
Features A — G are blue-shifted features and H — L red-shifted features, with respect to the systemic 
velocity, Vlsr = 10-8 km s“^. 

days, while others were detected during two or three epochs, possibly due to noisy observation 
at the third and fourth epoch. 

Figure ^ shows the obtained 3-D motions of the water-maser features. A clear separation 
between the blue-shifted and red-shifted maser features was found in the southeast and north¬ 
west directions, which indicate bipolar outflow. The former side of the outflow has a velocity of 
—8.68 km s“^ with respect to the systemic velocity, while the latter has a velocity of 5.45 km 
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Table 3. Obtained parameters of the water-maser features with proper motions. 



Position 


Relative velocity 



v^ 

drift 


Peak intensity 


Feature AX 

AY 

t-^X 



(7ij.Y 

V** 

Avz 

V 


Epoch 1 Epoch 2 Epoch 3 Epoch 4 


(mas) 

(mas) 

(km 


(km 

s-b 

(km 

s-b 

(km s" 

-Rr-b 


(Jy b 

sam ^) 


A ... 

48.23 

-68.58 

2.55 

2.35 

-4.00 

3.59 

-7.44 

0.73 

-1.37 

0.45 

7.87 

3.07 

1.43 


B ... 

46.20 

-64.82 

0.84 

7.45 

-6.34 

11.11 

-6.88 

0.42 

-1.27 

0.99 

5.48 

1.62 



c ... 

65.52 

-3.05 

7.66 

1.80 

-2.38 

3.06 

-5.41 

1.09 

-2.30 

0.45 

2.09 

6.32 

1.92 


D ... 

. -47.03 

36.59 

-0.25 

2.19 

0.93 

1.36 

-1.81 

0.76 

-0.68 

0.30 

1.45 

5.61 

2.32 

2.31 

E ... 

. -55.34 

42.34 

-4.64 

5.87 

4.11 

2.52 

-0.71 

0.99 

-1.55 

0.45 

3.64 

1.33 

1.47 


F ... 

53.42 

-37.67 

2.84 

2.07 

-6.2 

4.23 

-0.24 

1.58 

-4.32 

0.45 

2.62 

1.43 

0.86 


G ... 

. -62.24 

40.48 

3.68 

19.10 

4.78 

11.30 

0.42 

0.36 

-0.98 

0.99 

2.18 

0.73 



H ... 

. -67.80 

38.97 

-7.96 

7.02 

0.33 

1.30 

1.12 

0.95 

0.05 

0.30 

8.93 

3.91 

2.18 

2.97 

I ... 

. -72.40 

38.78 

0.67 

13.00 

-0.97 

5.08 

1.91 

0.73 

-2.52 

0.81 

2.19 

1.41 

0.64 


j... 

20.09 

46.81 

-0.62 

2.08 

2.95 

2.12 

2.72 

0.81 

-1.04 

0.45 

1.76 

1.29 

1.16 


Kt ... 

0.00 

0.00 

-5.54 

0.59 

4.46 

0.95 

4.23 

1.21 

-0.73 

0.30 

30.40 

30.07 

25.20 

20.35 

L ... 

. -0.46 

46.02 

0.77 

1.37 

2.32 

2.58 

6.56 

1.09 

-1.04 

0.45 

8.98 

2.17 

0.90 



Relative Doppler velocity with respect to the systemic velocity Vlsr = 10.80 km s 
^ Position reference feature at Vlsr = 15.03 km s~^ 


s“^. Features A, B, C, and F are within the blue-shifted side of the outflow, and spots G, H, 
and I are within the red-shifted side of the suggested bipolar outflow. Spots D and E possibly 
belong to the red-shifted side of the outflow, but have blue-shifted motions. 

3.1. Objective Analysis with VVCM 

To quantify the suggested bipolarity noted from hgure ^ we used the VVCM (Variance- 
covariance matrix) technique, which fully and objectively extracts the kinematic essentials, 
without assuming any particular model (|y00 [cite]cite.BloeBloemhof (2000)). The elements of 
the VVCM matrix are expressed as 
1 ^ _ 

( 1 ) 

n=l 

where the diagonal elements are the velocity dispersions. Here, i and j denote the three 
orthogonal space axes (R.A., Deck, and radial coordinates z), n is the n-th maser feature of 
the total N features, and the bar indicates the average over maser features. The VVCM matrix 
and the respective diagonalization for R Crt were obtained as follows (in units of km^ s“^): 


/ 18.55 

-7.35 

-8.10 ) 


( 35.82 

0 

0 \ 

-7.35 

15.97 

11.43 


0 

11.62 

0 

V -8.10 

11.43 

18.79 ) 


1 0 

0 

5.86 ) 


In hgure |^, the two bold vectors show two-dimensional projections of two unit eigenvec¬ 
tors with the largest eigenvalues of the VVCM diagonalization. The eigenvector V\ with the 
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Fig. 2. 3-D motions of maser features represented in scaled cones. The maser features are at the apex 
of the cones. The red and blue cones show red- and blue-shifted features, with respect to the systemic 
velocity. The inclination of the cone indicates the degree of the Doppler velocity with respect to the 
transverse velocity. The big plus-sign at the center is the model-calculated position of R Crt; the size 
of the plus-sign shows the positional uncertainties. The bold arrows are eigenvectors: Vi, the largest 
eigenvector (the VVCM outflow axis) and V 2 , the second-largest eigenvector. 

largest eigenvalue, Vi= 35.82 km^ s“^, is the major principal axis of the VVCM; the direction 
of the axis corresponds to the outflow axis. The axis of Vi lies at a P.A. of 136°, and is inclined 
at an angle of —39° with respect to the plane of the sky with the southeast lobe directed out 
of the page (towards the observer). The second-largest eigenvector V 2 , with the second-largest 
eigenvalue, ¥ 2 = 11.62 km^ s“^, is at P.A. of 68° east of north, and is inclined at an angle of 
24° with respect to the plane of the sky; the vector is directed into the plane. The collimation 
factor between the two largest eigenvalues, Vi and V 2 , is 3.08, and 1.98 between the smallest, 
V 2 and V3. On the other hand, a Monte Carlo simulation gave the following values: Vi= 70.6 
± 52.9 at P.A. 66°.5 ± 64°.9, an inclination of -4°.5 ± 23°.6, ¥ 2 = 24.2 ± 11.5 at P.A. 34°.2 
± 59°.3 and an inclination of 26°.8 ± 25°.4. This simulation gave the following values of the 
collimation factors: 3.1 ± 1.9 for ¥ 1 , ¥2 and 3.2 ± 1.8 for ¥ 2 , ¥ 3 . 
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We also analyzed the two-dimensional spatial variance-covariance matrix (SVCM). The 
SVCM matrix and diagonalization were obtained as follows: 

/ 2778.91 -1721.37 \ / 4130.36 0 

-1721.37 1937.82 / ^ \ 0 586.37 

Diagonalizing the matrix, one eigenvalue is larger than the other by a factor of 7, indi¬ 
cating a slightly elongated distribution. The axis of elongation lies at a P.A. of 128°, which is 
close to the value of 136° found for the velocity vector Vi. Despite the large uncertainties, the 
coincidence of these two model-independent axes supports the existence of a bipolar outflow of 
R Crt with a P.A. of 136° and an inclination of —39°. 

The smaller eigenvalues give a collimation factor of 1.98, the inequality of the eigenvalues 
suggests that the outflow of R Crt is azimuthally asymmetric around the outflow axis Vi. This 
means that transverse to Vi the cross section of the velocity dispersion ellipsoid perpendicular 
to the outflow axis Vi is elongated in the direction of eigenvector V 2 . From only the 12 maser 
features detected, it is difficult to state the signihcance of an azimuthal velocity asymmetry 
with an uncertainty of 98%. V 2 probably has a statistical bias that depends on the number of 
maser features. However, |e]cite.BloeBloemhof (2000) ([cite]cite.BloeBloemhof (2000)) found a 
collimation factor of 1.54 for water masers in W 49N, and suggested that magnetic helds may 
cause such azimuthal asymmetry. Due to flux freezing, ionized matter moves most easily along 
magnetic held lines and experiences magnetic drag in other directions. Therefore, how in a suf- 
hciently strong held may have an elongated cross section in velocity space, and V 2 may indicate 
the direction of elongation. |e]cite.Szy-aSzymczak et ah (1999) ([cite]cite.Szy-aSzymczak et al. 
(1999)) measured magnetic helds by observing linearly polarized OH 1667 MHz masers of R 
Crt in 1995. The results indicate that the magnetic held subtends a position angle of —30°, 
almost parallel to vector Vi. This may be a coincidence, since the magnetic held observation 
lasted for ~ 3 years before our observations. Therefore, we need simultaneous observations 
of a magnetic held and maser proper motions in order to conhrm the correlation between the 
magnetic helds and eigenvector Vi. 

3.2. Constructing the 3-D Kinematical Model 

In order to derive kinematic parameters of the CSE around R Crt, the 3-D motions 
of the features were htted to an expanding outhow model. First, the 3-D model assumes an 
expanding outhow where the velocity vectors depart radially from a common center. The details 
of the htting were described by ^]cite.lm-blmai et al.(2000) ([cite]cite.Im-bimai et al.(2000)). 
We estimated only the position and the systemic motion of the star as free parameters, but 
assumed the values of the distance to be d = 170 pc and the stellar velocity to be 10.8 km 
s“^. Table ^summarizes the estimated parameters in the model htting. Second, the expansion 
velocity was estimated by htting the data to the 2-D standard expansion model. This model 
traces an ellipse on a radius-W sr plot, and provides a reasonable description of the water maser 
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Table 4. Best-fit model for the maser velocity field for R Crt. 


Parameter Offset 


Velocity: 

Vqx (km s“^) . 5.4±3.9 

Voy*(km s"^) . -5.7±3.9 

Vbz^(km s“^) . 10.8 

Position: 

xq (mas) . —9±21 

yo* (mas) . 10±13 

RMS residual V . 0.47 


Relative value with respect to the position-reference maser feature. 

^ Assuming the systemic velocity: Vqz — 10.8 km s~^. 

CSE. In the case of the standard expanding ontflow model, the projected radius, r, of the maser 
feature at a velocity P — Vq is given by 

r = n[l-\(V-Va)IV,^,fYI\ (4) 

where Vq is the systemic velocity of the star, tq the CSE radius, and I4xp the expansion velocity 
at the CSE. Figure |] shows a radius- I4sr plot for R Crt and modeled ellipses. Least-squares hts 
of equation using the data (A — L features) gave average values of lAxp = 7.4 km s“^ and a 
CSE radius corresponding to tq = 2.6 x 10^^ m. The inner envelope hts to the parameters of r* 
= 1.3 X 10^^ m and I4xp = 4.3 km s“^. |e]cite.ColoColomer et ah (2000) ([cite]cite.ColoColomer 
et al. (2000)) estimated the radius of the water maser envelope of R Crt using the VLA as 2.1 
X 10^^ m with Vxp = 8 km s“^ at an assumed distance of 170 pc, which is consistent with our 
estimate with compact and bright features (c.f. ^00 [citejcite.Bo-dBowers, Johnston(1994)). 
The inner radius of the water-maser envelope emission is slightly larger than the typical radius 
of a dust shell with infrared emission (^00 [cite] cite.DancDanchi et al. (1994)). Thus, maser 
features are likely to he just above the dust shell where maser clouds are accelerated by radiative 
pressure (000 [cite]cite.DancDanchi et al. (1994) and references therein). The inner radius 
corresponds to a few stellar radii. The radius of the OH emission estimated using the MERLIN 
array (000 [cite]cite.Szy-aSzymczak et al. (1999)) is 3.6 x 10^^ m, where lAxp = 7.9 km s“^. 
|e]cite.Szy-aSzymczak et al. (1999) ([cite]cite.Szy-aSzymczak et al. (1999)) remarked that this 
OH envelope should be one of the smallest. Material in the CSE of R Crt is accelerated to reach 
a terminal velocity of R = 11.0 km s“^ (000 [cite] cite. ZuckZuckerman, Dyck (1986)), which was 
estimated from CO observations at a maximum radius of 3.1 x lO^'^ m (000 [cite]cite.Kaha- 
aKahane, Jura (1994)) at our assumed distance. Thus, the expansion velocity increases with 
distance. Note, however, that each of the different observations were performed at different 
epochs. 
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Fig. 3. Model fits assuming a standard expanding outflow. The radius of the outer edge of the water 
maser region is 100 mas, where the expansion velocity is 7.4 km s“^. The radius of the inner edge is 50 
mas, where the expansion velocity is 4.3 km s“^. 


The strength of the acceleration can be evaluated by comparing the observed expansion 
velocity, Kxp, with the escape velocity, 14 sc 5 given by 


hesc 


2GM, 


ro 


(5) 


where G is the gravitational constant and M* the mass of the star. Assuming a mass of R Crt 
of 1 Mq, at the inner radius (r* = 1.3 x 10^^ m, Resc = 14 km s“^, Kxp = 4.3 km s“^), 14xp 
is less than Rgsc; consequently, masers in the inner envelope are gravitationally bound to the 
star. At the outer radius (ro = 2.6 x lO^^m, Resc = 10 km s“^, Kxp = 7.4 km s“^), masers 
are marginally gravitationally bound. Analyses of the velocity helds of supergiants (e.g. VY 
CMa, S Per, and NML Cyg) suggest that material is, in general, bound to the star at the inner 
boundary, but unbound at the outer boundary (POO [cite]cite.YatesYates, Cohen (1994); ^00 
[cite]cite.RichRichards et ah (1996)). The strength of the acceleration in the case of R Crt may 
be less than those of supergiants. 

Finally, we discuss some relationships between the accelerations in the maser re¬ 
gions, the mass-loss rates of the central stars, and the asymmetry of the CSEs. The ac¬ 
celeration is parametrized by the logarithmic velocity gradient, e = d(ln n)/d(ln r) (|00 
[cite]cite.RichRichards et ah (1996)). For R Crt, within the outer limit of r = 100 mas and the 
inner limit of r = 50 mas, the logarithmic velocity gradient is e = 0.78. Table ^ summarizes 
the velocity gradients, the mass-loss rates, and the asymmetry of the CSEs in R Crt and in 
some AGB stars. Table shows that the asymmetry of CSEs seems to occur in stars with low 
velocity gradients (e < 0.8, ex., X Her, R Crt, NML Cyg, and VX Sgr), but not to occur in 
those with high velocity gradients (e > 1.2, ex., RT Vir and S Per). A large number of samples 
would clarify this correlation. 
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Table 5. AGB stars mass-loss rates and velocity gradients at the water-maser regions. 


Source 

name 

Type 

Mass-loss rate ^ 

10-® Mq yr-^ 

Velocity gradient ^ 

€ 

Morphology ^ 

X Her ... 

Semi-regular 

0.2 

__ (5) 

Bipolar 

R Crt ... 

Semi-regular 

1.0 

0.78 t®’ 

Bipolar 

RT Vir . . 

Semi-regular 

3.0 

3.50 

Spherical 

IK Tau . . 

Mira 

6.0 

0.84 (’■>* 

Spherical 

S Per .... 

Supergiant 

27.0 

1.20 

Spherical 

VY CMa 

Supergiant 

100.0 

0.55 

Ellipsoidal 

NML Cyg 

Supergiant 

180.0 

0.30 

Bipolar 

VX Sgr . . 

Supergiant 

270.0 

0.50 

Bipolar 


Upper limits. 

1 Mass-loss rates: ^jcite.Kaha-bKahane, Jura (1996) ([citejcite.Kaha-bKahane, Jura (1996)), ^jcite.Kaha- 

aKahane, Jura (1994) ([citejcite.Kaha-aKahane, Jura (1994)), 

Pejcite.YatesYates, Cohen (1994) ([cite]cite.YatesYates, Cohen (1994)) and references therein, 
^jcite.Eng-aEngels et al. (1983) ([citejcite.Eng-aEngels et al. (1983)). 

^ Velocity gradients: not available, this work, j^jcite.YatesYates, Cohen (1994) ([citeJcite.YatesYates, 

Cohen (1994)), 

l^jcite.RichRichards et al. (1996) ([citejcite.RichRichards et al. (1996)), 
^jcite.ChapChapman, Cohen (1986) ([citejcite.ChapChapman, Cohen (1986)). 

^ Morphology: Iniai H. private commnnications, ^ejcite.Bo-cBowers et al. (1993) ([citeJcite.Bo-cBowers 

et al. (1993)), 

^jcite.Rlch-bRichards et al. (1999) ([citejcite.Rich-bRichards et al. (1999)), yajcite.Rich- 
aRichards et al. (1998) ([citeJcite.Rich-aRichards et al. (1998)) 


3.3. Doppler Velocity Drifts 

One of the aims of this work is to measure the Doppler-velocity drifts of individual 
water-maser features. We found Doppler velocity drifts in the range from 14 = —4.32 km s“^ 
yr“^ (for feature F) to 14 = 0.05 km s“^ yr“^ (for feature H). Figure |T] and table ^ indicate 
that the blue-shifted maser features are apparently accelerated and the red-shifted decelerated. 
If the blue-shifted maser feature drifts imply true acceleration motions, the blue-shifted side 
of the outflow would be accelerated towards the terminal velocity, but the red-shifted features 
would not be accelerated. Analyses of the velocity gradients and the escape velocities suggest 
that water masers are accelerated in the envelope. It is expected that the red-shifted masers 
are decelerated if the driving source of the acceleration does not affect the red-shifted masers 
at the outer limit of the water-maser envelope. 

Remarkable velocity drifts were found for features F, the most negative-drifted, and 
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feature H, the most positive-drifted, F being blue-shifted, and H red-shifted with respect to 
the systemic velocity at 14sr = 10.8 km s“^. Since features F and H are near the systemic 
velocity, they move together with the star on the plane of the sky. Features F and H also have 
a relatively large proper motion on the celestial plane, F having 6.8 km s“^ and H having 8.0 km 
s“^ at the assumed distance. On the plane of the sky, features F and H are diagonally separated 
by r'U 100 mas, and the assumed stellar position lies almost at the center of the diagonal line 
between the two features. This line between features F and H coincidentally lies at P.A. ~ 
136°, which is almost parallel to the VVCM major axis. It is thought that different features 
do not switch on and off in a small area within the epochs of the observations, because of the 
simple structure of the maser features. We cannot reject the possibility that the drifts do not 
reflect the real motions of the circumstellar material. 

3.4- Distance to R Crt 

The distance to R Crt has been ambiguous, between 300 pc (000 [cite]cite.KnapKnapp 
et al.(1998)) and 170 pc (000 [cite]cite.Szy-bSzymczak, Le Squeren (1999)). For semi-regular 
variable stars, their intrinsic luminosity is unknown, and the distances can only be crudely 
estimated. Due to the small number of proper motions, no distance measurement based on 
statistical parallax and a model-htting method was available from our data. Instead, we made 
“spectrum-like prohles”, flux-density prohles against proper motions in the P.A. = 136° axis 
(direction of the VVCM eigenvector), and in the direction of the respective perpendicular axis 
P.A. = 46°. In hgure (b) represents the velocity prohles projected onto the directions in P.A. 
= 136°, and (c), that projected onto the P.A. = 46° axis. 

The spectrum coverages were 16.3 km s“^ in the radial direction, 14.8 km s“^ for P.A. 
= 136° and 10.8 km s“^ for P.A. = 46° axis directions at the assumed distance of 170 pc. 
It is likely that the expanding how is slightly blocked in the northeast—southwest direction. 
Thus, when assuming a distance of 170 pc, the spectrum-like prohles indicate that the velocity 
coverages in the radial and transverse directions are roughly equal. If a distance of 294 pc is 
assumed, the velocity coverages on the celestial plane would be double, and would imply an 
unreliable collimation of the bipolar outhow on the celestial plane. Therefore, a distance of 170 
pc is available for R Crt. Note, however, that the maser velocity coverage varies due to a time 
variation in the maser features. 

4. Conclusions 

We investigated the 3-D kinematics of water masers around the semi-regular variable 
star R Crt using data from multi-epoch observations with VLBI at short intervals. It has been 
verihed that most of these maser features existed for at least 80 days. We measured the three- 
dimensional velocity held, which indicates a possible bipolar outhow. A variance-covariance 
analysis applied to the velocity held supports the bipolarity of the outhow. The existence 
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Fig. 4. Spectrum profiles of water masers around R Crt. (a) Total-power spectrum obtained at the first 
epoch, j98090. Spectrum-like profiles derived from proper-motion measurements, (b) on the P.A. = 136° 
axis and (c) P.A. = 46° axis directions. The dashed lines at the bottom are the velocity widths, and the 
solid vertical line at the center is the systemic velocity. 

of bipolar outflow suggests that a Mira variable had already formed a bipolar outflow. We 
estimated the inner and outer radii of the water-maser region, which are consistent with those 
of other observations of R Crt. Water-maser emission in R Crt comes from regions where dust 
grains are condensed, and where material within the envelope is accelerated until reaching the 
terminal velocity. It seems that within the water-maser envelopes around R Crt, material is 
bound to the star in the water-maser regions, and marginally bound outside. We measured the 
Doppler-velocity drifts of the maser features. The blue-shifted maser features are apparently 
accelerated, while the red-shifted features are the opposite. If the accelerations of the blue- 
shifted masers are real, maser features are accelerating toward the terminal velocity. Making 
spectrum-like prohles from proper motions allowed us to verify the distance of 170 pc to R Crt. 
Similar observations are encouraged to directly measure the accelerations in CSE to investigate 
envelopes around stars and to clarify the asymmetric shapes of the CSEs of evolved stars. 
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